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O , Abstract. A new method for the determination of open cluster membership based on a cumulative effect is 

■ proposed. In the field of a plate the relative x and y coordinate positions of each star with respect to all the other 
stars are added. The procedure is carried out for two epochs ti and t2 separately, then one sum is subtracted from 
another. For a field star the differences in its relative coordinate positions of two epochs will be accumulated. 
For a cluster star, on the contrary, the changes in relative positions of cluster members at ti and t2 will be very 

1—1 1 small. On the histogram of sums the cluster stars will gather to the left of the diagram, while the field stars will 

■ form a tail to the right. The procedure allows us to efficiently discriminate one group from another. The greater 
the distance between ti and t2 and the more cluster stars present, the greater is the effect. The accumulation 

Q , method does not require reference stars, determination of centroids and modelling the distribution of field stars, 

^ • necessary in traditional methods. By the proposed method 240 open clusters have been processed, including stars 

ly-j ' up to m < 13. The membership probabilities have been calculated and compared to those obtained by the most 

, commonly used Vasilevskis-Sanders method. The similarity of the results acquired the two different approaches 

is satisfactory for the majority of clusters. 

^ . Key words, open cluster - membership probabilities- proper motion 
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m . 

\0 1. Introduction 

vn 

Various methods based on the analysis of positions, proper motions, radial velocities, magnitudes and their combina- 
I tions have been proposed to determine the members of open clusters, 
f — The first mathematically rigorous procedure for determination of open cluster membership was developed by 
' Sanders (|197ip with a statistical analysis of proper motions. It is also the most widely used method. Sanders's approach 
^ [ is based on the model of overlapping distributions of field and cluster stars in the neighborhood and within the region 
•""j ■ of visible grouping of stars, introduced by Vasilevskis et al. (1958). 

r> I Vasilevskis's model implies that proper motion dispersion of cluster members is caused by observational and 
' measurement errors assumed to be normally distributed. Thus the distribution of cluster stars is represented by a 
bivariate normal frequency function. The dispersion of field stars is due to not only the errors referred to above, but 
also to peculiar motion and differential galactic rotation. Therefore the field star distribution is not expected to be 
random, but rather to have a preferential direction and not normal distribution. However, in a first approximation 
a bivariate normal ellipsoidal distribution function was assumed for field stars, the major axis of the ellipse being 
parallel to the galactic plane. 

Thus, Sanders's equations contain 8 unknown parameters to be determined: the number of cluster members, x 
and y components of two centroids, one circular and two elliptical dispersions. This system is solved by a maximum 
likelihood method (Sanders ll97ip . The probability of stars being cluster members is calculated by frequency functions 
with determined parameters. 

Slovak (1977) tested the Vasilevskis-Sanders method by modelling the proper motion distribution in the surround- 
ings of a cluster. He proved the uniqueness and convergence of solutions of the system, provided that errors are 
represented by a Gaussian distribution and open clusters have no significant internal motion. But if there is noticeable 
motion within a cluster or it rotates then the above method fails. 

Cabrera-Cafio and Alfaro (|1985p improved the numerical techniques for obtaining the above parameters. McNamara 
and Schneeberger p978p showed that the final probabilities could be influenced by various weight groups. Zhao and 
He (|1990p provided a method for treating data with different accuracies. 
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However all these improvements did not treat the problem of star distribution model on which the membership 
probabilities are based. Even if the hypothesis of a normal distribution of field stars is realistic for some clusters, the 
centroids of field and cluster stars sometimes are too close to be well discriminated. The parametric model also fails 
when the cluster member-to-field star ratio is small. The Vasilevskis- Sanders method does not work in the case of 
significant internal motion in a cluster or its rotation. 

To overcome some of the problems arising from the parametric Vasilevskis-Sanders method, especially the star 
distribution modelling, Cabrera-Caho and Alfaro l|1990p developed a more general, non-parametric method of mem- 
bership determination. Here no assumptions were made about the nature of cluster and field star distributions and 
allowed for the use of photometric data. Each cluster needed careful individual study. 

So, the actual distribution of cluster and field stars, especially the latter, may not be fitted by a Vaselevskis-Sanders 
model; the centroids for the two groups may be too close to be distinguished; the reliability of the method depends on 
the cluster-to-field star ratio; in the case of significant internal motions or rotation of a cluster the traditional method 
fails. 

In the next section we introduce a new method for the discrimination of cluster stars from surrounding field 
stars based on a cumulative effect using positions and proper motions. We do not try to overcome all the difficulties 
of traditional methods, though we offer another approach that enlarges the statistical distance between the two 
populations - cluster members and field stars - by revealing the group of stars with the least relative velocities. The 
advantage of this method is that no assumption is made about the distribution of field stars and determination of 
centroids is avoided. However, in order to determine probabilities we still have to assume a normal bivariate distribution 
for clusters. Another advantage of the method is that reference stars are not necessary: the discrimination of cluster 
members is most effective by rectangular coordinates. These features of the method allow us to increase the statistical 
distance between the two populations. The most noticeable advantage of the accumulation method is its ability to 
reveal dynamic structures within the clusters if there are any. 



2. The method 

We consider a cluster as a gravitationally bound system with a bulk motion relative to surrounding galactic field stars. 
Field stars are characterized by more or less random velocities. Therefore the cluster members must show regular 
motion, while the field stars move realtively irregularly. 

We consider rectangular coordinates of stars in the field of a plate (or CCD array), measured in two epochs ti and 
^2. The relative x and y positions of each star with respect to all the other stars in the area are added. If the two 
epochs are well separated, then for a field star the differences of the relative coordinate positions in the two epochs 
will be accumulated. For cluster stars, the changes of relative positions of physical members at ti and ^2 must be 
very small and their sum will be minimal. The more cluster stars, the bigger is the effect. If one draws a histogram 
of differences, the stars with minimal relative velocities (which we consider more likely to be cluster members) will 
gather to the left of the diagram, while the field stars will form a tail to the right. Thus, the procedure allows us to 
efficiently discriminate one group from another, avoiding the use of reference stars. 

For each star the following quantities are calculated: 



N 

i=i 

N 

i=i 



N 



where iV is a total number of stars (on the plate or CCD array). The primed quantities belong to the first epoch, the 
double primed to the second one. 

Then the differences of S'-,.. and S'J.. are calculated: 

The next step is to draw the histograms for 6Sxi ^^vi separately. According to our assumption, cluster stars 
will concentrate within the first bin with minimal {SSxi , SSy- ) . 

In order to define the probabilities of belonging stars to a cluster, the first one or two most populated bins are 
chosen. Then the whole procedure is repeated for the stars from the selected bins using the formulae (1) and (2), but 
dropping the moduli this time, that is: 



N 



N 

K = E(^^ " y'^ 



(4) 
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N N 

Again, the histograms are drawn. To determine the probabihties of stars being physical members, we use the normal 
ellipsoidal distribution function as the best fit for a cluster in the first approximation. So, the obtained diagrams are 
fitted with Gaussian curve by a least square method: 

/(a;) = ^=e ^ "x \ f{y) = ' "v ' , (6) 

where A is an area of the normal distribution curve, cr is a double dispersion. 
The probabilities by x and y coordinates for each star are calculated as: 

= e-'( — > , py^e ' , (7) 

The resulting probability is: 

P^Px*Py (8) 

If there is no "raw material" (i.e. rectangular coordinates) available, but there are proper motions instead, then 
(1), (2), (4) and (5) arc modified as follows: 

JV N 
N N 

The rest of the procedure remains the same. In this article we apply this modified version of the method to the 
clusters from Bias's pOOll I2002p catalogues and compare the results. 

As will be shown later, the accumulation method can reveal various dynamic groups in a cluster, if there are any. 



3. Results 

The accumulation method was applied to 240 open clusters of the Bias catalogues, including the stars up to m < 13. 
The clusters in the Bias list were processed by a Vasilevskis-Sanders approach. In order to compare our results, we 
introduce a degree of similarity for each cluster as: 



where are probabilities by Bias, - ours, n - number of stars in a cluster. For C a histogram was drawn (see 
Fig- Hi- The smaller the C the closer are the results obtained by the two different methods. 

As one can see from the histogram, the similarity of probabilities is quite satisfactory for the great majority of 
clusters: 162 out of 240 have C < 25, 65 clusters have C in the range between 25 and 45 and 13 of them 45 < C < 70. 
The latter group have the most serious discrepancies - almost an anticorrelation. 

We choose two cases to show: Collinder 121 (Fig. [2]), with C — 7, total number of stars n = 179 - good agreement, 
and NGC 3228 (FigE]), with C = 50, n = 32 - disagreement. 

In Fig. [2] and FiglHl the histograms in the upper row are for \Sfj,a\ and \5^i3\ according to the formulae (9), those 
in the middle row - for 6fia and Sfj^p by the formulae (10). In the figure on the third row the individual probabilities 
by Bias (denoted by circles) and by the accumulation method (asterisks) are presented. 

The probabilities for the cluster Collinder 121 (Fig. [5]) reveal close similarity, as for 2/3 of clusters, while those of the 
cluster NGC 3228 (Figl3]) differ considerably from Bias's results. In the latter case the anticorrelation of probabilities 
is striking: the stars we claim as cluster members are field stars by Bias and vice versa. In this case (number of stars 
32) the probabilities are on average smaller than ours. Our probabilities are somewhat more pronounced: cluster stars 
have high probabilities, while those by Bias are low. According to Bias, it may not be considered as a cluster at all, 
since only two stars have membership probabilities greater then 50(only 16 out of 32 have V and B — V quantities 
measured), as shown in Fig. ID 13 stars out of 16 lie along one sequence confirming that this is an open cluster with 
high probability. In the table below we give V and B — V magnitudes and probabilities by Bias (Pd) and by the 
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accumulation method {Pa) for each star. Star No 17 is claimed to be a cluster star by Dias, while we claim both stars 
Nos 17 and 18, which are far from the main sequence of the diagram, to be field stars. 

The reasons for the disagreements between our results in this case, as well as in those few cases with poor accordance 
could be several: the low cluster-to-field star ratio, or the presence of more then one dynamic structure or large 
observational errors. 

Pleades (Figl5])is an example of dynamic structure or a rotation. The histograms for Smuai (upper diagram) and 
SmuSi (lower diagram) are shown, the quantities calculated using (9) and (10). Two peaks on both diagrams are clearly 
seen, which reflects either two dynamic structures or rotation of the cluster. 

4. Conclusion 

In the presented method of accumulation we assume that the cluster members moving in the Galaxy as a whole 
have similar velocities, while field stars show a wide range of velocities. The procedure - adding close velocities while 
compensating for disperced ones - should enhance the assembling of cluster stars, while distributing more sparsely the 
field ones, due to the cumulative effect. Thus this method effectively enlarges the statistical distance between physical 
members and field stars, so that a membership as well as a non-membership is more pronounced. This effect is better 
for larger cluster-to-field star ratios. For field stars no particular distribution is assumed, thus the centroids are not 
determined. However, when calculating probabilities, for physical members a normal bivariate distribution function 
is assumed. For rectangular coordinates no reference stars are needed in the procedure. The most interesting feature 
of the accumulation mathod is its ability to reveal more then one group of velocities, which has been shown with the 
example of Pleades. 

The modified accumulation method was applied to 240 clusters from Bias's list. The probabilities calculated by 
the accumulation method showed satisfactory agreement with those obtained by the Vasilevskis-Sanders method for 
the majority of clusters. The poor agreements or disagreements can be ascribed to low cluster-to-fleld star ratios, or 
multiple dynamic structures. 

The accumulation method can be further expanded using other variables like radial velocities and photometric 
quantities. 
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Star No 


V 


Bv 


Pd 


Pa 


2 


8.14 


0.96 


10 


84 


3 


9.09 


-0.1 


12 


72 


4 


8.43 


-0.09 


11 


84 


5 


8.19 


-0.02 


9 


81 


6 


7.90 


-0.03 


16 


57 


8 


9.37 


-0.05 


10 


92 


9 


9.34 


0.14 


27 


43 


11 


9.97 


0.05 


19 


49 


12 


9.31 


-0.06 


12 


85 


13 


9.03 


-0.06 


9 


90 


15 


11.26 


0.33 


8 


98 


16 


10.96 


0.13 


35 


25 


17 


9.54 


1.99 


69 


17 


18 


9.79 


1.56 








19 


7.97 


-0.09 


15 


57 


20 


10.17 


0.11 


9 


92 
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Fig. 1. Histogram for C - the degree of similarity of the results obtained by Dias and the accumulation method (11). 
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Fig. 2. CoUinder 121: the first row - the histograms for |(5/iQ,| and \Sfifj\ (formulae (9)); the second row - histograms 
for 5^a and 5^(3 (formulae (10)); the third row - comparison of individual probabilities by Dias (circles) and by the 
accumulation method (asterisks) 
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Fig. 3. NGC3228: the first row - the histograms for \6^a\ and ISn^] (formulae (9)); the second row - histograms for 
Sfia and SiJ,j3 (formulae (10)); the third row - comparison of individual probabilities by Dias (circles) and by the 
accumulation method (asterisks) 
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11.5 




B-V 



Star V B-V P_D P_A 
No 



2 


8,14 


0,96 


10 


84 


3 


9,09 


-0,10 


12 


72 


4 


8,43 


-0,09 


11 


84 


5 


8,19 


-0,02 


9 


81 


6 


7,90 


-0,03 


16 


57 


8 


9,37 


-0,05 


10 


92 


9 


9,34 


0,14 


27 


43 


11 


9,97 


0.05 


19 


49 


12 


9,31 


-0.06 


12 


85 


13 


9,03 


-0.06 


9 


90 


15 


11.26 


0.33 


8 


98 


16 


10.69 


0.13 


35 


25 


17 


9,54 


1.99 


69 


17 


18 


9,79 


1.56 








19 
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-0.09 


15 


57 


20 


10.17 


0.11 


9 


92 



Fig. 4. The color-magnitude diagram for NGC3228 
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Fig. 5. The Pleades (Mellotte22): histograms for djion (left) and Sji^i (right). 



